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Abstract. We give a leisurely introduction into mathematical diffrac- 
tion theory with a focus on pure point diffraction. In particular, we 
discuss various characterisations of pure point diffraction and common 
models arising from cut and project schemes. We finish with a list of 
open problems. 

1. Introduction 

Quasicrystals were discovered by their unusual diffraction properties [55J . 
Subsequently, quite some mathematical effort has been devoted to diffraction 
theory of aperiodically ordered models (see corresponding parts in [101 ED 
[5*4} [61] and references therein). In fact, there is substantially more work 
done than could sensibly be reviewed here. For this reason we will focus 
below on pure point diffraction and model sets. These topics seem to be 
both particularly relevant and conceptually fairly well understood. We note 
in passing that model sets were introduced by Meyer [40] quite before the 
dawn of quasicrystals. They became a standard in the physical literature. 

Recent results give a diffraction theory for measures on locally compact 
abelian groups [6] [Tl [TT] [38] . Here, however, we will restrict our attention 
to point sets in Euclidean space in order to keep the exposition as simple 
as possible. The much studied topic of mixed spectra and random systems, 
see e.g. [3 [SJ [131 SSI EQ] and references therein, will only appear in the last 
section on problems. 

For earlier survey type articles on mathematical diffraction we refer to 
|28[ [12] . A detailed introduction into mathematical diffraction theory is 
given in the lecture notes [2]. Another somewhat introductory text from the 
point of view of dynamical systems can be found in [37J. 

The article is organised as follows: In Section [2] we introduce the frame- 
work of mathematical diffraction theory and define and discuss its key quan- 
tity, the diffraction measure. Section [3] is concerned with characterisation of 
pure point diffraction. In Section [3] we discuss cut and project schemes and 
model sets. Finally, in the last section we discuss some open problems. 



2. The framework of mathematical diffraction theory 

Diffraction theory for crystals has a long history. Discussions can be found 
in many textbooks, e.g. [18]. The first systematic treatment of mathematical 
diffraction theory for aperiodic order is due to Hof [24l [25] . Here, we discuss 
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the basic setup. The mathematical key quantity turns out to be a measure, 
the so called diffraction measure. It models the outcome of a diffraction 
experiment. 

We will model the solid in question by a subset A of Euclidean space R . 
We will assume that A is distributed in a regular way with points being not 
too close and not too far. More precisely, we will assume that A is Delone, 
i.e. 

• there exists an r > such that different points of A have distance at 
least 2r (" A is uniformly discrete" ) , 

• there exists an R > such that no point of Euclidean space has 
distance bigger than R from A (" A is relatively dense"). 

Thus, clearly, A has infinitely many points. In order to understand the 
diffraction setup for the infinite A it is helpful to first consider the case of a 
finite set F of scatterers. In this case the intensity is the function given by 
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exp(-i(x - y)£). 

x,y<=F 



The analogous expression for the infinite A reads 

x,y£A 

This sum is heavily divergent, as there are infinitely many terms of modulus 
1. This is not a mathematical problem only. There is a physical reason 
behind this divergence: The intensity of the infinite solid is infinite. The 
correct object to be considered is the normalised intensity per unit volume. 
It is defined as 

(!) T A = lim -TB-rlAnBn- 

n-*oo \r> n \ 

Here, B n denotes the ball around the origin with radius n and the modulus 
denotes the volume. Note that A n B n is finite by our assumption on A and, 
hence, lAnB„ is a well defined function. There are issues in (pQ) in which sense 
the limit is meant (it is the vague sense) and whether it exists at all (we 
may have to pass to a subsequence). However, we will skip these somewhat 
technical details. 

The quantity I a is known as diffraction measure. It describes the outcome 
of a diffraction experiment. The values £ with Ia({£}) > are called Bragg 
peaks. The value of Ia({€}) is called intensity of the Bragg peak. 
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Now, that we have defined I a we may wonder how to calculate it. To this 
end let us expand the definition of I a- This gives 

1 . ,.1 
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Here, <5 X denotes the unit point mass at x and F is the Fourier transform. 
This calculation shows that I a is the Fourier transform of 

1 1 x,ye4n.B„ 

This object ja is known as autocorrelation or Patterson function (though it 
is a measure). 

In a theoretical sense the investigation of diffraction is thus reduced to 
the following two step procedure: (1) Calculate the Patterson function. (2) 
Take its Fourier transform. 

While it is possible to carry out this procedure for various examples, it is 
far from obvious how to do this in the general case. For this reason a lot of 
effort has been put into finding general criteria for the investigation of point 
spectrum. In this context, the basic questions are the following: 

• When is I a a pure point measure? 

• Where are the Bragg peaks? 

• What are the intensities of the Bragg peaks? 

Some answers to these questions will be discussed in subsequent sections. 

3. Pure point diffraction 

On the conceptual level two approaches to pure point diffraction have 
been developed in the last ten or so years. These are given, respectively, by 

• considering the associated dynamical systems, or 

• investigating almost periodicity properties. 

In this section we will be concerned with these two approaches. 

3.1. Pure point diffraction via dynamical systems. When dealing 
with disordered systems in statistical mechanics it is quite standard not 
to consider a single object but a whole ensemble of objects exhibiting the 
"same type of disorder". In the range of aperiodic order the same reason- 
ing can be applied. It suggests to consider not a single set A but rather 
the ensemble of all subsets r of Euclidean space which have the same local 
structure as A. This ensemble will be denoted by fi. It can be made into 
a compact topological space [HI EHJ [29] but we will not worry about this 
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here. Instead we will note the following crucial feature of 17: If T belongs 
to 17 then so does its translate by t G M> N given as 

t + r = {t + x ■. x g r}. 

The obvious reason is that T and t + T have the same local structure. Thus, 
17 together with the translations gives a dynamical system. We will write 
(17, R ) to denote this structure. 

In many situations 17 comes with a canonical translation invariant mea- 
sure m. In these cases it is possible to provide an autocorrelation 7 = j m 
by a closed formula invoking only m and not using a limit. This was first 
realized by Gouere |22j . While his further considerations use the theory of 
stochastic processes and Palm measures, the formula for j m can rather di- 
rectly be given. Following Baake/Lenz [6] we obtain for the application of 
the measure 7 m to the continuous compactly supported function tp 




where a is a function on 1 N with f o~(t)dt = 1. In this way the autocorrela- 
tion is defined without a limit. It can then be shown that it actually agrees 
with the corresponding limit almost surely [221 [6] . 

Whenever 17 is equipped with an invariant measure m it is further possible 
to consider the associated space L 2 (l7,m) of square integrable functions on 
17. The translations on 17 induce a unitary representation of M. N on L 2 (17, m) 
via 

(T t f)(r) = f(t + r) 

for t G R . A function / G L 2 (17, m) is called an eigenfunction to the 
eigenvalue £ if 

T t f = exp(itO/ 

for all t G R . If there exists a basis of L 2 (£2, m) consisting of eigenfunctions, 
the dynamical system (Q,M. N ) is said to have pure point spectrum. The 
main result on pure point diffraction in the context of dynamical systems 
now reads as follows. 

Theorem 1. Let (17, R-^) together with an invariant measure m be given. 
Then, the following assertions are equivalent: 

(i) The Fourier transform I m of 7 m is a pure point measure. 

(ii) The dynamical system (17, R-^) with the measure m has pure point 
spectrum. 

In this case the group generated by the Bragg peaks ({£ : I m ({€}) > 0}) is 
the group of eigenvalues of the dynamical system. 

This theorem gives a characterisation of pure point diffraction in terms 
of the dynamical system. Moreover, it provides further information on the 
position of the Bragg peaks. The implication (ii) ==>■ (i) of the theorem 
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established in [19] has been a basic tool in establishing pure point diffrac- 
tion for concrete classes of models |24[ [56| l58~l [53] . The converse implication 

(1) => (ii) allows one to set up a perturbation theory for pure point diffrac- 
tion in the context of dynamical systems. This has been carried out in [7] 
(see [17] for related material as well). A particular application is the study 
of deformed model sets, as done earlier in [141 123] . 

The theorem is the outcome of cumulated efforts over many years: The 
analogue result for one dimensional subshifts was proven by Queffelec in [49] . 
For the dynamical systems we consider here, the implication (ii) ==> (i) is 
due to Dworkin [19] with later modifications by Schlottmann [56] and Hof 
|25j (see [20] for related material as well). The full equivalence was proven 
for systems satisfying a certain regularity assumption known as finite local 
complexity by Lee/Moody /Solomyak in [33 j . For the systems we consider 
(and even more general ones) the full equivalence was then shown by Gouere 
[23] , It is possible to leave the class of point processes and to go to measures 
instead as discussed in Baake/Lenz [6] and then also in Lenz/Strungaru [39] . 
The statement on the eigenvalues is implicit in the proof of [33]. An explicit 
formulation can be found in [6]. 

The theorem rises the question whether the spectrum of the dynamical 
system is always given by the diffraction spectrum. This turns out to be 
wrong as was shown by van Enter/Mi§kisz [20J. More precisely, the dy- 
namical spectrum always contains the diffraction spectrum [19]. However, 
it maybe strictly richer than the diffraction spectrum as shown by examples 

The theorem does not answer the question on the intensities of the Bragg 
peaks. The basic physical intuition concerning this issue is that the intensity 
of the Bragg peak at £ should be given by 

(2) I A ({0) = lim 

n— >cx 

Validity of this equation is sometimes discussed under the header of Bombieri/Taylor 
conjecture. In fact Bombieri/Taylor use validity of this equation in their 
work [151 \W\ without any further justification. By now validity has been es- 
tablished by direct arguments for primitive substitution systems by Gahler/Klitzing 
[21] and for regular models sets |24[ 156]. In fact, these results follow from 
a conceptual approach to ([2]) via continuity of eigenfunctions and uniform 
Wiener /Wintner type results. This has been developed by Lenz [36 . The 
conceptual approach itself is hinted at in |25[ [28] . The necessary continu- 
ity of eigenfunctions is proven for model sets in [251 [56] and for primitive 
substitutions by Solomyak in [59J. The corresponding Wiener /Wintner type 
results had already been studied by Robinson [52] . 

3.2. Pure point diffraction and almost periodicity. In some way or 
other notions of almost periodicity have been around in the study of pure 
point diffraction for quite a while. In particular, there is work of Solomyak 
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[57j providing a connection and a discussion of Lagarias [28] asking for con- 
nections. There is even a characterisation of pure point diffraction in terms 
of almost periodicity in the work of Queffelec in [49] for the (somewhat dif- 
ferent) situation of symbolic dynamics. Furthermore, Meyer's work on what 
is now known as Meyer sets [40] and subsequent discussions, see e.g. [42] . 
have a very almost periodic flavour. Still it seems fair to say that only with 
the rather recent work of Baake/Moody [llj and Gouere [23] the strength 
of this connection became apparent. 

A continuous function / on R N is called almost-periodic if for any e > 
the set of its e-almost-periods 

{teR N : \\f(--t)-f\\ 00 <e} 

is relatively dense in H. . Here, || • ||oo denotes the supremum norm. Similar 
notions for measures exist and give the concept of strongly-almost-periodic 
measure and norm-almost-periodic measure. While we do not want to con- 
cern the reader with the technical definitions here, we would like to note 
that norm almost periodicity is substantially stronger than strong almost 
periodicity. 

It is also possible to develop a concept of almost periodicity for sets. This 
is done under the name of Bohr/Besicovitch almost periodicity by Gouere 
[23j (see corresponding questions of Lagarias in [28J as well). 

Theorem 2. Let A be a Delone set whose autocorrelation ja exists. Then, 
the following assertions are equivalent: 

(i) The Fourier transform I a of ja is a pure point measure. 

(ii) ja is a strongly- almost-periodic measure. 

(iii) A is Bohr/Besicovich almost periodic. 

For sets A whose autocorrelation is supported on a uniformly discrete set, 
the result is proven by Baake/Moody [llj . In the general form stated above 
this theorem is due to Gouere [23] . 

As will be discussed in the next section, one is particularly interested 
in the case of Meyer sets A. These are Delone sets with the property that 
A— A is uniformly discrete. For such sets it makes sense to define an e-almost 
period or better a statistical e-almost period t G R N with 

j(A \(A + t)U(A + t)\ A)) n B n 
hmsup — — ; < e. 

n— >oo \-*-*n\ 

Here, ft denotes cardinality. For Meyer sets Baake/Moody [11] have the 
following result. 

Theorem 3. Let A be Meyer with autocorrelation ja- Then, the following 
assertions are equivalent: 

(i) The Fourier transform La of ja is a pure point measure. 

(ii) ja is a norm- almost-periodic measure. 
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(iii) For any e > the set of statistical almost-e-periods of A is relatively 
dense. 

Let us emphasise that the actual setting and results of [11] are consider- 
ably more general than discussed in this theorem. In particular, [11] deals 
with weighted point sets on locally compact abelian groups. Moreover, it 
gives a natural connection between pure point diffraction and cut and project 
scheme (see below for details). 

4. Cut and project schemes and model sets 

The two most prominent classes of mathematical models for aperiodic 
order are models coming from primitive substitutions and models coming 
from cut and project schemes. The latter are often discussed under the name 
of model sets (see e.g. |42} 143] for further discussion and references). They 
provide also standard examples discussed in the physical literature. There 
is a wealth of results on model sets and cut and project schemes. Here, we 
focus on the following issues: 

• explicit computation of I a for regular model sets, carried out by Hof 
[24~1 [25] (see later generalisations [561 CD] as well) , 

• existence of a lot of point diffraction for general sets associated to 
cut and project schemes, shown by Strungaru [60] . 

• a natural connection between cut and project schemes and pure point 
diffraction, discovered by Baake/Moody [11], and then further ex- 
plored in [32j H3 [8] , 

• characterisation of primitive substitutional sets with pure point diffrac- 
tion as model sets, due to Lee [30] (see the work of Barge/Kwapisz 
[?. 126] for an analogous one dimensional result in a slightly different 
context). 

The first result justifies mathematically the corresponding parts of the 
physical literature. The second result shows that order in the sense of Meyer 
condition implies order in the sense of a large point component in the diffrac- 
tion spectrum. The third result (or rather the corresponding circle of ideas) 
shows that cut and project sets arise very naturally within the framework 
of pure point diffraction for Meyer type sets. The final result shows that the 
"other class of examples" viz primitive substitutions is not really a different 
class when it comes to models with pure point diffraction. 

Let us now start by shortly recalling the framework of a cut and project 
scheme. Besides the physical space 1R a cut and project scheme has two 
further ingredients. These are a further space and a lattice. The further 
space is known as internal space, perpendicular space or reciprocal space. It 
will be denoted by H. It does not need to be an Euclidean space. It suffices 
if it is a locally compact Abelian group. The lattice is denoted by L. It is 
a lattice in K. x H. Its projections to the physical space and the internal 
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space will be denoted by L and L* respectively. A precise definition of cut 
and project scheme now runs as follows. 

A cut and project scheme over ~K N consists of a locally compact Abelian 
group H, and a lattice L in R N x H such that the canonical projection 
7r : R N x H — ► M. is one-to-one between L and L := vr(L) and the image 
7r int(-C') of the canonical projection 7Ti nt : R N x H — > H is dense. Given 
these properties of the projections ir and 7r int , one can define the *-map 
(.)*: L — ► H via x* := (7r int o (ttIl)" 1 ) (x), where (tt\l)^ 1 (x) = ir^ 1 (x) HL, 
for all x G L. We summarise the features of a cut and project scheme in the 
following diagram: 

u u u 



L ► L* 

We will assume that the Haar measures on M. N and on H are chosen in such 
a way that a fundamental domain of L has measure 1 . 

Given a cut and project scheme, we can associate to any W C H, called 
the window or atomic surface, the set 

X (W) := {x G L : x* G W} 

A set of the form t + A (W) is called model set if the window W is 
relatively compact with nonempty interior. 

Theorem 4. Let A be a Delone set. Then, the following assertions are 
equivalent: 

(i) A — A is uniformly discrete. 

(ii) A is a subset of a model set. 

(hi) There exists a finite set F with A — A C A + F. 

The equivalence of (i) and (ii) is due to Meyer [3D] and Moody |42| . The 
equivalence of (i) and (hi) is due to Lagarias [27]. The sets characterised in 
the previous theorem are known as Meyer sets. Note that all three conditions 
appearing in the theorem can be understood as indicating long range order 
in form of a weak lattice type condition. Various further characterisations 
can be found in the literature |42j . 

4.1. An explicit formula for Let a cut and project scheme (M. N , H, L) 
be given. Let also a sufficiently nice window W in H be given. Sufficiently 
nice means roughly speaking that the window is not a fractal. More precisely, 
we require W to be compact with non empty interior and and a boundary 
of measure zero. In this case, one can calculate explicitely the diffraction 
measure I a [24^ [56] . We need the dual lattice L 1 - of L given by 

L x := {(k,u) G x H : e m u(l*) = 1 for all (1,1*) G L}. 
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Let L° be the set of all k £ M for which there exists u € H with (k, u) £ L . 
This set is sometimes known as reciprocal lattice. It can be shown that there 



If H happens to be an Euclidean space as well, the formula for reads 

I f w e*p( ik *y)dy\ 2 - 

4.2. A lot of point diffraction for Meyer sets. In this subsection we 
highlight the following result of Strungaru |60| . 

Theorem 5. Let A be Meyer with autocorrelation 74. Then, j2 has a 
relative dense set of Bragg peaks. 

As mentioned already this result can be understood as linking two notions 
of long range order. The result is rather general as it does not assume any 
further regularity properties of the point set in question. 

4.3. A natural cut and project scheme. It is a fundamental insight of 
Baake/Moody that any set with a sufficiently almost periodic autocor- 
relation comes with a natural cut and project scheme. The required almost 
periodicity of the autocorrelation in turn is equivalent to pure point diffrac- 
tion whenever the autocorrelation is supported on a uniformly discrete set. 
This ties pure point diffraction and cut and project schemes (within the 
context of Meyer type sets) . 

A crucial step in the argument of [11] is to use the autocorrelation function 
to introduce a topology on point sets. Taking completions with respect to 
this topology then yields the internal space. In this way, cut and project 
schemes lie at the crossroads of two topologies: the local topology and the 
topology coming from the autocorrelation function. This point of view is 
further developed in [321 l4"4"l l4"6] . 

These results allow one to derive the characterisation of pure point diffrac- 
tiveness for Meyer sets given in Theorem [3] above. They can also be used 
to characterise the dynamical systems arising from regular model sets. This 
has been discussed by Baake/Lenz/Moody [8]. As [IT] provides a cut and 
project scheme, the main task in [8] is to construct and study the window 
using properties of the dynamical system. The basic connection between 
the dynamical system and the cut and project scheme is given by the so 
called torus parametrisation [4j [56] . The torus parametrisation turns out to 
be strongly determined by properties of the eigenfunctions. Two properties 
of the eigenfunctions are central. These are continuity of eigenfunctions and 
their separation properties. 



exists a unique group homomorphism * : L° 



H such that 
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As a byproduct one obtains a characterisation of lattices within Meyer 
sets with pure point diffraction. The work [8] also plays an important role 
in the investigation of substitution systems discussed next. 

4.4. Substitutions with pure point diffraction. As mentioned already 
the most important classes of examples for aperiodic order are model sets 
and sets arising from primitive substitution. By their very construction mod- 
els arising from primitive substitutions have a strong form of self-similarity. 

Here, we will shortly discuss a remarkable result of Lee |30] (see [?] as 
well) relating primitive substitution sets to model sets within the context 
of pure point diffraction. In some sense the result shows that one can not 
get away from model sets when dealing with pure point diffraction. More 
precisely, Theorem 5.5 of [30] gives in particular the following. 

Theorem 6. Let A be induced by a primitive substitution with finite local 
complexity. Then, the following assertions are equivalent: 

(i) A has pure point dynamical spectrum. 

(ii) A is an inter model set. 

Let us point out that the notion of a set associated to a primitive sub- 
stitution requires some care. One possibility is to consider tilings arising as 
fixed points of a primitive substitution. Then a point is chosen in each tile 
in a consistent way (i.e. so that points for tiles of the same type are in the 
same relative position in the tile). Thinking of the points in different classes 
of tiles as marked with different colours we obtain a coloured set. Such a 
set is underlying the statement of the previous theorem. We refrain from 
giving precise definitions for a coloured set but continue to further explain 
the statement of the theorem. 

The condition (i) is equivalent to a (suitably defined) notion of pure point 
diffraction for coloured sets. The condition (ii) also has to be understood for 
coloured sets. By definition an inter model set agrees with a model set up to 
points induced by the boundary of the window. Instead of using tilings one 
can directly deal with substitutions for point sets, see Theorem 5.3 in |30| . 
In this case the additional assumption of legality of certain clusters has to 
be imposed. 

The proof of these results winds together different strings of reasoning. 
One such string concerns so called coincidence conditions. They give criteria 
for a primitive substitution to have pure point diffraction. Building up on 
earlier work [31|. [34] Lee gives a new coincidence condition allowing for a 
characterisation of pure point diffraction for primitive substitutions. This 
characterisation is in fact part of the main result. It is done under the 
assumption that the set in question is Meyer. The second string then is a 
result of Lee/Solomyak [35] showing that a primitive substitution with pure 
point diffraction must be Meyer. As a third ingredient Lee uses the recent 
results of Baake/Lenz/Moody [8] providing a characterisation of regular 
model sets in terms of the associated dynamical systems. 
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5. Open questions 
In this section we present various issues and questions for further research. 

• Geometric implications of pure point diffraction. 

The phenomenon of pure point diffraction seems to be fairly well under- 
stood within the context of Meyer sets. This poses the question whether 
pure point diffraction in itself together with mild geometric restrictions (as 
e.g. finite local complexity and repetitivity) forces the Meyer property. For 
primitive substitutions this has been answered affirmative by Lee/Solomyak 
[35] . The general case seems to be open. A particular instance of this type 
of issue is the question which properties single out the lattices within the 
sets with pure point diffraction. For further discussion of these and related 
issues we refer to the article of Lagarias [28] . 

A further issue in this context is the validity of Bombieri/Taylor conjec- 
ture (discussed above) for larger classes of examples. 

Another question concerns entropy. Of course, entropy should somehow 
vanish for models with aperiodic order. Indeed, sets with pure point diffrac- 
tion and further regularity can be shown to have vanishing topological en- 
tropy [9]. On the other hand, there are natural examples of sets with pure 
point diffraction exhibiting positive topological entropy. Such an example is 
given by the set of visible points as shown by Pleasants [IB]. It seems that 
this is related to cut and project schemes with a window with a "thick" 
boundary. 

• Mixed spectra and random systems. 

The understanding of mixed spectra is very much at the beginning. Let 
us illustrate this by considering two extreme cases: On the one hand there 
are primitive substitutions models. These models carry a lot of order by 
their very construction. Still, not all primitive substitutions have point 
spectrum, let alone pure point spectrum. Thus, mixed spectra go well with 
a very rigid order structure. This is a conceptual issue in the understanding 
of order as encoded by spectral properties. On the other hand there are 
random systems. Random systems based on lattices exhibit a tendency to 
have mixed spectra with a pure point component due to the lattice and an 
absolutely continuous component due to the randomness. While this is well 
confirmed in examples [13] a general treatment is not available yet. 

Actually, random systems and substitutions are not that far apart in 
terms of diffraction. More precisely, as discussed by Hoeffe/Baake [5j it is 
possible to have a primitive substitution system with the same diffraction 
as a random system. 

One reason that diffuse spectra are not as well understood as point spectra 
is that there does not seem to be a good dynamical interpretation. 

• Homometry and inverse problem. 

The above considerations have been concerned with the direct problem 
i.e. to determine the diffraction given the solid. Of course, the real problem 
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is the inverse problem. In mathematical terms this amounts to describing all 
configurations leading to the same diffraction. This is known as homometry 
problem. In this context one may ask for properties shared by all solutions 
to the inverse problem as well as for further restrictions making the solution 
unique. 
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